Bright tripartite entanglement in triply concurrent parametric oscillation 
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Entanglement is a central property of quantum me- 
chanics. In particular, continuous-variable (CV) entan- 
glement is at the core of the Einstein-Podolsky-Rosen 
(EPR) paradox and is also important for quantum 
teleportation, quantum information, and quantum cryp- 
tography 0, U- Bipartite entanglement is now read- 
ily producible experimentally, and there has been some 
progress in the production of tripartite entangled beams, 
where entanglement was obtained by mixing squeezed 
vacua with linear optical elements In our proposal, 

strong entanglement is created in the nonlinear interac- 
tion itself and, as we will demonstrate, is present above 
the operating threshold, where the entangled outputs are 
macroscopically bright beams. We confirm that the out- 
put fields satisfy criteria for measurable CV tripartite 
entanglement [g and derive two types of experimental 
EPR [l| criteria which are applicable to this system. We 
then show analytically and numerically that the proposed 
system demonstrates these properties, both above and 
below threshold. Note that the tripartite CV entangled 
state thus created tends towards a GHZ state in the limit 
of infinite squeezing, but is analogous to a W state for 
finite squeezing [g. 

It is worthwhile to discuss the relation of the multipar- 
tite CV entangler studied in this paper to the multipar- 
tite CV entangler proposed by van Loock and Braunstein 
and realized recently || . The latter is comprised of N 
squeezers (some of which may be replaced by vacuum in- 
puts if non-maximal entanglement is acceptable) and an 
V-input-port interferometer. The experiment therefore 
requires stabilization of the frequencies and output pow- 
ers of N optical parametric oscillators (OPO), as well as 
of the N optical paths of the interferometer. In contrast, 
our system is based on a single OPO pumped by a fre- 
quency comb (which can be easily provided by a single 
CW laser, frequency- or amplitude- modulated, or by a 
single, stable, pulsed laser) and there is no interferome- 
ter at the output. However, the nonlinear medium inside 
the OPO must provide what we have named concurrent 



nonlinear interactions, i.e. simultaneously phase-matched 
three-field mixings where each field produced comes from 
two of the nonlinear processes. (Fig. H). 




(c) (d) (e) 

FIG. 1: Graphical representation of multipartite CV entan- 
glers for N = 5. Each dot represents a cavity field mode and 
the lines represent their respective nonlinear couplings. The 
solid lines in a graph represent, indifferently, all down- or all 
up-converting nonlinear interactions; the dashed lines then 
represent the interactions opposite to the solid line ones, i.e. 
all up- or down-converting, respectively, (a): ideal concur- 
rence, yields maximal (GHZ) entanglement in the limit of in- 
finite squeezing, (b): van Loock and Braunstein's ideal case: 
yields maximal (GHZ) entanglement in the limit of infinite 
squeezing, (c): nonideal concurrence due to an incomplete 
set of interactions, (d): nonideal concurrence due to an unde- 
sirable degenerate interaction of the same kind as the nonde- 
generate ones, (e): nonideal concurrence due to a "leaking" 
interaction. In (c,d,e), entanglement is limited to a number 
of modes much smaller than N. 

For example, the following Hamiltonian 

H = ih (xJ>i2o[at + X^2344) + H c -i W 
where the bij denote the pump fields, describes a concur- 



2 



rence, even though the multipartite entanglement avail- 
able from this interaction is not very strong [^, 0. If 
a\ and a 2 are swapped in the second term 0, 0, 0] , 
the available tripartite entanglement increases a little, 
although it is still not as strong as for the system we 
analyse in the present work. In the graphical represen- 
tation of Fig. the case of Eq. (JJJ corresponds to three 
vertices connected by only two edges, i.e. an incomplete 
triangle, analog to Fig. [He). 

In the general case, the following concurrence generates 
a set of N entangled modes ai,...,iv @ : 



N 

H(N) = ifyc'yi'y] hja\a] + H.c. 

i=l j>i 



(2) 



where we have taken all nonlinear coupling constants 
equal to x f° r simplicity. The pump fields bij were 
assumed to be classical undepleted fields in Ref. H, 
whereas they will be quantized in the present analysis. 
The analysis of Refs. shows that the Hamiltonian 

Hrw) leads to perfect multipartite entanglement and cor- 
responds to Fig. nja). However, one must: (i) have all 
possible pair couplings, (ii) avoid all couplings outside 
the set of N entangled modes, and (Hi) carefully sup- 
press all degenerate interactions (i.e. terms of the form 
bi+kA-k <v) since just a few of these destroy entangle- 
ment very quickly. Examples of such undesirable con- 
figurations are given in Fig.^c,d,e). Their quantitative 
effects are presently under study and the results will be 
published elsewhere |l4|. 

It is important to note that degenerate interactions do 
not destroy entanglement if they have the opposite effect 
from the nondegenerate ones, i.e. upconverting if the non- 
degenerate ones are downconverting, or vice versa. This 
is expressed in the Hamiltonian by opposite signs for the 
nondegenerate and degenerate terms @ and, in Fig. ^ 
by solid and dashed lines, respectively. In this case, if all 
possible degenerate interactions are present with equal 
coupling strength [Fig. Gib)], then the situation becomes 
exactly that of the van Loock and Braunstein proposal 
||. This is, however, clearly impossible to implement in 
a single OPO and all degenerate interactions must there- 
fore be suppressed in our proposal. 

In the absence of any degenerate interaction, the rela- 
tive signs between nondegenerate Hamiltonian terms do 
not necessarily decrease the amount of entanglement, 
even though they do change the particular entangled 
state that is created. For example, for N=3, The Hamil- 
tonian _ff( 3 ) (J2J, with classical and equal pumps, is 



H, 



(3) 



ihx 



H.c. 



(3) 



and admits eigenstates of the following (unnormalized) 
form, in the amplitude quadrature basis (and to constant 



amplitude shifts left) 

1^(3)) = J \x)i\x) 2 \x) 3 dx 



(4) 



where p a is a constant parameter. When one sign is 
changed in Hrg\, for example, 



H[ 3) = ih x 



-a[a 2 + a 2 a 3 + a 3 a 1 



H.c. 



the eigenstate becomes of the form 



1^(3)) = / \p)i\ph\ -P>3 dp. 



(5) 



(6) 



which is a state orthogonal to ^(3) but is still maximally 
entangled. Another example is 



H (3) 



ihx 



-a\a,2 — S2Q.3 + Sgaj 



H.c. 



(3)/ 



I - x) 2 \x)z dx. 



(7) 
(8) 



It is thus easy to see that no entanglement is lost 
with sign changes for nondegenerate interactions between 
three modes. It is reasonable to assume it might still be 
the case for N > 3 but we haven't yet investigated this 
more general situation. 

Finally, an important fundamental point is that, as 
recently shown by Braunstein [l^, there exists a formal 
mathematical connection between our entangler and that 
of van Loock and Braunstein: a linear algebraic trans- 
formation called the Bloch-Messiah reduction projects 
our system onto theirs. (In technical terms, the Bloch- 
Messiah reduction is a particular case of singular-value 
decomposition of the Bogoliubov transformation matrix 
of the system.) 

The ideal concurrences of Fig.^a) would be extremely 
difficult to obtain using birefringently phase-matched 
nonlinear optics; however, the advent of quasi-phasc- 
matching has radically changed this situation, as ex- 
plained in Ref. , and an experimental demonstration of 
three simultaneous nonlinear interactions, involving the 
same set of wavelengths in the same nonlinear crystal, 
has recently been achieved 01 . A concrete and natu- 
ral way of implementing Ht N \ experimentally is to have 
the entangled modes' labels i,j correspond to the eigen- 
frequencies of an optical resonator, such as the resonant 
cavity of an OPO. In that case, additional mode quantum 
numbers are provided by the optical polarization, a two- 
dimensional eigenspace. It is then possible to experimen- 
tally design quasi-phase-matched concurrences for closed 
sets of modes, up to N — 4, using setups such as the one 
described below. (See Refs. 0, f° r the experimental 
details.) These setups are free of any additional terms 
deleterious to entanglement. In order to achieve sc aling 
of entanglement to a larger number N > 4 of modes [l5j , 
it appears mandatory to use a higher-dimensional quan- 
tum number than polarization. The angular momentum 
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of light created in Laguerre-Gaussian beams might 
be an interesting possibility to explore. 

We now describe the experimental system analyzed in 
detail in this article. A schematic of our system is given 
in Fig. showing the three quasi phase-matched pump 
inputs, which interact with the crystal to produce three 
output beams at frequencies ujq, lo\ and oj 2 , with the in- 
teractions selected to couple distinct polarisations. Mode 
1 is pumped at frequency and polarisation (loq + wi, y) to 
produce modes 4 (u>q, z) and 5 (oj\, y), mode 2 is pumped 
at (u>i +u>2,y) to produce modes 5 and 6 (uj 2 ,z), while 
mode 3 is pumped at (2u>i, z) to produce modes 6 and 4. 
Note that x is the axis of propagation within the crystal. 
It is important to note that no other interaction involving 



the usual route 01 we obtain the master equation 
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FIG. 2: Schematic of the experimental setup. Pump lasers 
drive three modes (denoted by dots, squares and circles), with 
suitable frequencies and polarisations, which are down con- 
verted to three other output modes by the crystal. Note that 
the physical beams will overlap in the crystal, the separation 
being intended here for clarity. 

the three output modes is phase-matched in this system. 
To briefly confirm the feasibility of this scheme: first, 
it is easy to fulfill the cavity resonance conditions for 
such a closed set of 3 coupled modes, the problem being 
similar to the well-known doubly resonant type-II OPO. 
Second, type-I cluster instabilities, which are expected 
above threshold for the zzz interaction, can be avoided 
by making the frequencies different enough for dispersion 
to be significant in the nonlinear medium, or by adding a 
mode filter inside the cavity. We will now turn to the full 
quantum analysis of the concurrence-based OPO, below 
and above threshold. 

The Hamiltonian for the six-mode system is 



Htot = Hi IC 



Hh 



(9) 



where the rotating frame interaction Hamiltonian is 

flint = ^(xiaiaJ4 + X2a2a^4 + X3a 3 44) + H - c -' ( 10 ) 

with the Xj representing the effective nonlinear couplings, 
and the free and damping Hamiltonians have their usual 
forms. The external pump fields (but not the intracav- 
ity ones) are treated as classical. Note that the analy- 
sis of Ref. ||, in contrast, was limited to H lnt with all 
pump modes being classical and no cavity whatsoever. 
Its conclusions were therefore applicable only to the op- 
tical parametric amplifier case, not to the OPO above 
threshold, one of the situations treated here. Following 



dp 
at 



Hint, P] 



(11) 



where the Lindblad superoperator T>j\p] = 2ajpa J 

OjCijp— pa'jaj is obtained by tracing over the density ma- 
trix for the usual zero-temperature Markovian reservoirs. 
We assume that all intracavity modes are resonant with 
the cavity, although detuning can easily be added. We 
will treat all of the high (low) frequency modes as having 
the same cavity loss rate 7 (k). In what follows we set 
7i = 72 = 73 = 7) an d 74 = 75 = 76 = noting that 
the down conversion efficiencies are not necessarily 
identical 

We now map the master equation onto a Fokker-Planck 
equation for the positive-P function Note that we 

must use the doubled phase space of the positive-P in 
order to ensure positive semi-definite diffusion. We can 
then establish a correspondence between stochastic am- 
plitudes ctj (and ctj~) and mode operators dj (and at) 
respectively, ctj and ci^ being independent complex vari- 
ables. The full quantum dynamics of the system can now 
be found by solving the Ito stochastic differential equa- 
tions 



{-jai + Ex - XiO-iCi^) dt, 
(-7a 2 + E 2 - X2Ct 5 a6) dt, 
(-70:3 + E 3 - X3aAa 6 ) dt, 
(—Ka 4 + xiaia^ + X3a3<a;g ) dt 
+ VxIaT dWi (t) + Vxi^i dW 3 (t) , 
(— K015 + xi&iot^ + xiuiol^) dt 
+VX^ dW 2 (t) + VxI^T alW{ (t) , 

(-K06 + X2d2c4 + X3"3a4 ) dt 

+Vx^ dw;(t) + vx^i dw;(t), (12) 

and also the equations obtained by interchange of ctj with 
ct - , and dWj(t) with dWj + z{t). The six independent 
Gaussian complex noises are completely determined by 
their non-vanishing correlator 



da 1 
da 2 
da 3 
da\ 

da 5 

dae. 



dW*(t)dWj(t') = §ij§{t - t') dt. 



(13) 



Defining the vector a = (a\, af, a 2 , a 2 , . . . , ae, o.q) t , 
we determine the stability of the system by linearising 
the equations about their semi-classical steady state so- 
lutions (denoted by a). We first drop the noise terms in 



Eq. El so that a 



and solve for the steady state 



solutions. The linearised equations for the fluctuations 
about the steady state 5a = a, — a. then form a mul- 
tivariate Ornstein-Uhlenbeck process [2(|. The resulting 
matrix equation is 



dSa = A8adt + BdW, 



(14) 
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where dW is a vector of independent real noises, B is 
the noise matrix of Eq. ^] with the steady-state values 
inserted, and 



Ai -A 2 
(A* 2 ) T A 3 



(15) 



where A\ = —7/6 



and 
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, (16) 



(17) 



For simplicity in what follows, we will assume that all 
the nonlinearities are equal, i.e. Xj = X- This can be 
achieved via quasi-phase-matched interactions in period- 
ically poled ferroelectrics || , and means that we may set 
the effective coupling rates equal by equating the pump 
field amplitudes, Ej = E. 

If we treat all the pump fields as real, all the signal 
fields are also real since the fields must satisfy the phase 
condition </>4 + </>5 = </>4 + (/)$ = 05 + 0g = 0. As with the 
standard optical parametric oscillator there is an oscilla- 
tion threshold, which for this system occurs at the pump 
field strength 



= 7 «/2 X , 
below which the stable solutions are 



"(1,2,3) = E/j, 



*(4,5,6) 



= 0. 



(18) 



(19) 



At the threshold the fluctuations are not damped and 
therefore a linear fluctuation analysis is not valid. If 
the pumping is increased further, modes 4, 5 and 6 be- 
come macroscopically occupied, while the high frequency 
modes saturate at their threshold value. There is also a 
relationship between the driving fields that must be satis- 
fied to have a steady state, which in the case of equal non- 
linearities, means that all the pump amplitudes must be 
equal and hence all the signal amplitudes are equal |2lj . 
The resulting steady state solutions 



"(1,2,3) 



k/2x, a (4 ,5,6) - ^{E-E^)/x, (20) 



enable us to investigate the quantum correlations of the 
system by studying fluctuations around the steady state, 




FIG. 3: Spectra for the entanglement criteria of Eq. 1221 for 
equal pump amplitudes and nonlinearities in each mode. In 
all figures \ = 0.01, « = 1, 7 = 10 and the pump values are 
E = 0.9E th (solid line) and E = l.LE th (dash-dotted line). 
All quantities plotted in this article are dimensionless. The 
dashed lines represent the limiting zero-frequency threshold 
value of 2/9 and the tripartite correlation value of 4, which is 
obviously violated. 

from which we will also find the measurable extra-cavity 
fluctuation spectra poj . 

We will first describe the measurable quantities which 
may be used to experimentally verify that this system ex- 
hibits true multipartite entanglement. We define quadra- 
ture operators for each mode as 



Xj 



Y, 



-i(aj 



(21) 



so that [Xj, Yj] = 2i and the Heisenberg uncertainty prin- 
ciple requires V(Xj)V(Yj) > 1, which sets the vacuum 
noise level. Conditions which are sufficient to demon- 
strate bipartite entanglement for Gaussian variables are 
well known [2^ ■ These have been generalised to tripartite 
entanglement by Van Loock and Furusawa 0, without 
making any assumptions about Gaussian statistics. Us- 
ing our quadrature definitions, these conditions give a set 
of 3 inequalities 



V{Xi - Xj^) + V(Y 4 +% + %) > 4, 



(22) 



where i, j € {4, 5, 6} and V(A) = (A 2 ) - (A) 2 . The simul- 
taneous violation of any two of these conditions proves 
genuine tripartite entanglement for the system. An un- 
depleted pump analysis of the interaction Hamiltonian 
(Eq. I1U|) reveals that these are the maximally squeezed 
quadrature combinations, so that we can expect optimal 
entanglement for these variables @. 

The Einstein, Podolsky and Rosen (EPR) paradox 
stems from their 1935 paper [l|, which showed that local 
realism is not consistent with quantum mechanical com- 
pleteness. Optical quadrature phase amplitudes have the 
same mathematical properties as the position and mo- 
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FIG. 4: Output fluctuation spectra of the inferred quadra- 
tures for three mode EPR correlations, from measurements 
on one or two modes. Below threshold (E — 0.9E th , solid 
line) the spectra coincide. Above threshold (E = l.lE t , 
dash-dotted line) two mode inference gives better violation 
than one mode inference, but the difference is indistinguish- 
able on this scale. 

mentum originally considered by EPR and, as demon- 
strated by Reid [2j|, this property can be used to define 
inferred variances which, when entangled sufficiently, al- 
low for an inferred violation of the uncertainty princi- 
ple. This demonstrates the EPR paradox and was exper- 
imentally studied by Ou et al., who found clear agree- 
ment with quantum theory [24j |. Following the approach 
of Ref. [23, we assume that a measurement of the Xj 
quadrature, for example, will allow us to infer, with some 
error, the value of the X k quadrature, and similarly for 
the Yj^k quadratures. Minimising this error, we can de- 
fine variances of the inferred quadratures, the products of 
which demonstrate the EPR paradox when they appear 
to violate the Heisenberg uncertainty principle. 

Although tripartite entanglement is present in the sys- 
tem, it does not appear that a demonstration of the EPR 
paradox is possible via the standard two mode approach. 
There are two ways to demonstrate paradoxes of this na- 
ture, both involving information about all three modes. 
One example is to use quadrature measurements on one 
mode to infer values for the joint state of the other two. 
This is mathematically equivalent to inferring, for ex- 
ample, the center of mass position and momenta of two 
particles by measurements on a third. In this case we 
find, using mode i to infer properties of the combined 
mode j + k, for example, the inferred variances 



V mf (X J ±X k ) = V(X 3 ±X k ) 



V{X i ,X j ±X k ) 



V{Xi 



(23) 

with similar expressions holding for the Y quadratures. 



In the above V(A,B) = (AB) - (A)(B) and the i,j,k 
can be any of 4, 5, 6 as long as they are all different. 
The product of the actual variances always satisfies a 
Heisenberg inequality so that there is an experimental 
demonstration of a three mode form of the EPR paradox 
whenever 



V mf {X J ± X k )V mJ (Yj ± Y k ) > 4 



inf/ 



(24) 



is violated. Another option is to make measurements on 
the combined mode j + k to infer properties of mode i. 
This leads to the inferred variances 



V mf (X i ) = V(X i )- 



V(X u X 3 ±X k ) 



V(Xj ± X k ) 



and a demonstration of the paradox whenever 
V in f (Xi)V inf (YA > 1. 



(25) 



(26) 



is violated. 

The experimentally accessible quantities are the nor- 
mally ordered fluctuation spectra corresponding to the 
inequalities defined in Eqs. [221 El [23 f° r which the 
same relationships hold. These are defined as Fourier 
transforms of the time-normally ordered operator covari- 
ances, 

Sfa&Mu) = ?T(: (ii(t)-<ii(*)»(4(0)-(4-(0))) :), 

(27) 

which are related to the measurable output spectra, 
S out [Ai, Bj](ui), using the standard input-output rela- 
tionships for optical cavities 25] ■ In the results presented 
here, we will treat the output mirror as the only source 
of damping, which is a common approximation in theo- 
retical quantum optics. For the inequalities of Eq. 11' I'D 
we define the fluctuation spectra 

7<f H = S out [Xi - X^} + S out [Y A + %+ %], (28) 

and use the abbreviation S out [Aj\ = S^Aj, Aj]. Sim- 
ilarly, for the variances in Eqs. 1241 1261 we obtain vari- 
ances of fluctation spectra (S m '(u})) which must violate 
the same inequalities in the frequency domain. 

We can find relatively simple analytic expressions for 
the correlations of Eq. EHlin the case where all the pump 
amplitudes and nonlinearities are equal. We denote the 
sub-threshold solutions (Eq. I19fl by a = (5(1,2,3) j an d the 
super-threshold solutions (Eq. I2TJI) by (3 = a^ 5 ^ 6 y The 
symmetry means that J|^*(£j) = 7|^*(o;) = /^(w) = 
1^(00), where — (+) denotes the spectra below (above) 
the oscillation threshold. The measurable spectra then 
take the form 
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7° ut H = 5 



8xna[7(xa) 2 + lOx^a + A(lu 2 + k 2 )] 
{(xa + k) 2 + uj 2 } {(2xa + k) 2 + U) 2 } ' 

Ak 2 {lo 2 + 7 2 )[76(x/?) 4 + (x/3) 2 (100 7 k - 56w 2 ) + (lo 2 + 7 2 )(16w 2 + 43k 2 )] 



(4(x/3) 2 + 2 7 k - lo 2 ) 2 + {2k + 7 ) 2 lu 2 (2(x/3) 2 + 3 7 k - 2lo 2 ) 2 + (3k + 2 7 ) 2 uj 2 



(29) 
(30) 



The solutions for these quantities are shown for pump 
field amplitudes of 0.9E th and l.lE th in Fig. Below 
threshold the vacuum outputs exhibit near total violation 
of the van Loock-Furasawa entanglement criteria near 
zero frequency and approach the uncorrelated limit of 5 
for large frequency. Above threshold the outputs form 
bright beams and maintain a large violation of the in- 
equalities near zero frequency, which implies strong tri- 
partite entanglement. We note that at E = l.lE th the in- 
tensities of the three low- frequency modes are twice those 
of the high frequency modes. Fig.0]shows strong inferred 
violations of Eas. 1241 1^1 for the same degree of entangle- 
ment shown in Fig. |3 indicating that above threshold 
the system can be used to demonstrate three mode EPR 
correlations with bright output beams. 

In conclusion, we have shown that concurrent intra- 
cavity x nonlinearities can be used to produce strongly 
tripartite entangled outputs, both above and below the 
oscillation threshold. The above threshold entanglement 
is found for macroscopic intensities, demonstrating a 
bright tripartite entanglement resource. We have also 
shown how this device can be used to perform multi- 
mode demonstrations of the EPR paradox. We expect 
that this proposal will be easier to stabilise than devices 
which rely on separate OPOs and multi-port interferom- 
eters. 
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